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We report a direct lattice QCD calculation of the strange nucleon electromagnetic form factors 
Ge and G%i in the kinematic range 0 < < 1.2 GeV^. For the first time, both Gf and G%i are 

shown to be nonzero with high significance. This work uses closer-to-physical lattice parameters than 
previous calculations, and achieves an unprecedented statistical precision by implementing a recently 
proposed variance reduction technique called hierarchical probing. We perform model-independent 
fits of the form factor shapes using the ^-expansion and determine the strange electric and magnetic 
radii and magnetic moment. We compare our results to parity-violating electron-proton scattering 
data and to other theoretical studies. 


The nucleon electromagnetic form factors describe how 
electric charge and current are distributed inside protons 
and neutrons, and are therefore among the most impor¬ 
tant observables characterizing these building blocks of 
ordinary matter. Because nucleons contain only up and 
down valence quarks, these two quark flavors dominate 
the electromagnetic form factors. Isolating the small con¬ 
tributions from the other quark flavors is a significant 
challenge for both experiment and theory, but is of fun¬ 
damental importance for our understanding of the struc¬ 
ture of protons and neutrons, and of the nonperturbative 
dynamics of QCD. After the np and down quarks, strange 
quarks are expected to give the next-largest contribution 
to the electromagnetic form factors. The cross section of 
elastic electron-proton scattering used to extract the form 
factors is dominated by photon exchange, which probes 
the sum of all quark-flavor contributions weighted ac¬ 
cording to their electric charges. However, by analyzing 
the small parity-violating effects arising from interference 
with Z-boson exchange, the strange-quark contribution 
to the electromagnetic form factors can be isolated [1, 2]. 
The available experimental results, which focus on mo¬ 
mentum transfers in the vicinity of 0.2 GeV^, are con¬ 
sistent with zero but constrain the relative contribution 
of the strange quarks to be within a few percent [3-14] . 

Ab-initio calculations of the nucleon electromagnetic 
form factors and G\j of an individual quark flavor 
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FIG. 1. Nonperturbative quark propagator contractions ob¬ 
tained from Eq. (1) by performing the path integral over the 
quark fields: connected (left) and disconnected (right). 


q (see, e.g., Ref. [15] for the definitions) are possible us¬ 
ing lattice QGD. The form factors can be extracted from 
Euclidean three-point functions of the form 

y- e-m'-(-y)g-.p.(y-x) (Np{z)V^^{y)N^{x)) , (1) 

z. y 

where N is an interpolating field with the quantum num¬ 
bers of the nucleon, E/ = is the vector current for 
quark flavor q, and p, p' are the spatial momenta of the 
initial and final states. In the three-point function (1), 
performing the path integral over the quark fields leaves 
a path integral over the gauge fields, which contains the 
product of the fermion determinants and the nonpertur¬ 
bative quark propagator contractions illustrated in Fig. 1. 
The connected contraction arises only for q = u,d and 
is numerically large, while the disconnected contraction 
is present for all quark flavors (and is the origin of the 
strange-quark contribution to the electromagnetic form 
factors). The disconnected quark loop in Fig. 1 has the 
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form 

T!^ = -J2 Tr y)] , (2) 

y 

where Dq is the lattice Dirac operator, and the trace is 
over color and spin indices. The numerical computation 
of the propagator D~^{y, y) for all spatial lattice points 
y using standard methods is prohibitively expensive for 
lattices of realistic size, and therefore most lattice calcu¬ 
lations of nucleon form factors have been restricted to the 
connected contractions (and therefore q = u,d) only. In 
Ref. [16], the strange-quark contribution was estimated 
by combining experimental data for the complete elec¬ 
tromagnetic factors with lattice QCD results for the con¬ 
nected u- and d-quark contractions. This method relies 
on a delicate cancellation between large quantities, and 
is therefore limited in its statistical precision and rather 
susceptible to systematic errors in the lattice calculation. 

A direct lattice QCD calculation of the strange-quark 
contribution can be performed by evaluating the discon¬ 
nected loop in Eq. (2) stochastically. An unbiased esti¬ 
mate is given by 

1 ^ 

E E ^Uy)rMy), (3) 

n=l y 

where are suitable noise vectors (for example with ran¬ 
dom components G Z 2 ) with support on the time slice 
?/o, and ipn are the corresponding solutions of the lat¬ 
tice Dirac equation, Dqtpn = ^n- Previous unquenched 
lattice QCD calculations of the strange nucleon electro¬ 
magnetic form factors using variants of this approach 
can be found in Refs. [17] and [18]. The calculation of 
Ref. [17] was performed on a 16^ x 32 lattice with dimen¬ 
sions (1.9fm)^ X (3.9 fm), and at rather heavy u, d quark 
masses corresponding to pion masses in the range 600-840 
MeV. For these parameters, the strange magnetic form 
factor was found to be negative with a significance 
of 2-3(7 in the region ^ 1 GeV^, while the strange 
electric form factor G|; was found to be consistent with 
zero [17]. The authors of Ref. [18] used an anisotropic 
24^ X 64 lattice with dimensions (2.6 fm)^ x (2.3 fm) and 
a pion mass of 416 MeV. The corresponding results for 
G^ and have large uncertainties and are consistent 
with zero. 

In this Letter, we report a new direct lattice QCD cal¬ 
culation of and G% with closer-to-physical parame¬ 
ters and with very high statistical precision. For the first 
time, both G^ and Gg are shown to be nonzero with 
high significance and in a wide range of Q^. We also 
obtain similarly precise results for the analogous discon¬ 
nected light (equal for up and down) quark contributions, 
which need to be included in precision lattice calculations 
of the total electromagnetic form factors of the proton, 
and are therefore of relevance for the proton charge ra¬ 
dius puzzle [19, 20]. Our calculation is performed on a 





FIG. 2. The first 16 Hadamard vectors for hierarchical prob¬ 
ing on a 4^ lattice; red = -fl, black = — 1. 


large 32^ x 96 lattice of dimensions (3.6 fm)^ x (10.9 fm) 
and includes 2-1-1 flavors of dynamical sea quarks, im¬ 
plemented using a clover-improved Wilson action. The 
up and down quark mass corresponds to a pion mass of 
317 MeV, and the strange-quark mass is consistent with 
the physical value (determined using the “rjs” rnass [21]) 
within 5 percent. The unprecedented statistical precision 
is achieved as follows: (i) we use 1028 gauge-field config¬ 
urations and compute the three-point function in Eq. (1) 
for 96 different source locations, x, on each configuration, 
and (ii), we use a novel variance reduction method, hier¬ 
archical probing [22], to evaluate the disconnected quark 
loops T^. Similarly to dilution [23, 24], this method is 
based on the observation that the fluctuations in Eq. (3) 
due to the random noise vectors originate from the off- 
diagonal elements of D~^{x,y), which decay with the 
Euclidean distance \x — y\. By partitioning the noise vec¬ 
tors into multiple vectors with support only on subsets 
of the lattice sites, the variance can be reduced (at the 
cost of additional solutions of the lattice Dirac equation). 
While dilution is based on a fixed partitioning scheme, 
hierarchical probing allows one to continuously increase 
the level of partitioning, eliminating the variance in or¬ 
der of importance while reusing the results from prior 
levels [22]. This is achieved using a special sequence of 
Hadamard vectors, which have values ±1 on the lat¬ 
tice sites. Examples of are shown in Fig. 2. In Eq. (3), 
we make the replacement 


Zn 0 (4) 

where ^ is a single noise vector and © denotes the 
element-wise Hadamard product [22] . As mentioned ear¬ 
lier, we use noise vectors with support only on selected 
time slices, and we therefore perform three-dimensional 
Hierarchical probing. We use V = 128 Hadamard vec¬ 
tors, which eliminates the variance from neighboring lat¬ 
tice sites up to distance 4. For the electromagnetic 
form factors, we observe a variance reduction by approx¬ 
imately a factor of 10 compared to the traditional noise 
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method with the same N (at equal computational cost). 
Our calculation employs complex Z 2 x Z 2 noise and also 
uses color and spin dilution for the noise vectors 

We use the lattice conserved vector current [generaliz¬ 
ing Eqs. (2, 3) for currents that are not site-local] includ¬ 
ing an improvement term [25, 26] to reduce 0(a) lattice 
artifacts, where a is the lattice spacing, with the improve¬ 
ment coefficient set to its tree-level value. We indepen¬ 
dently vary the momenta p and p' to obtain 59 different 
values of Q^. For clarity in plots we bin the nearest points 
that have the same (p' — p)^ and Jp'^ — P^j, and thus dif¬ 
fer in by less than 0.006 GeV^. The three-point func¬ 
tion (1) receives contributions from the desired nucleon 
ground state and from excited states; the excited-state 
contributions decay exponentially faster with the source- 
sink separation, \zo — xq], than the ground-state contri¬ 
bution. We use \zo — xq] « 1.14 fm for our main results, 
and estimate the remaining excited-state contamination 
using a second separation of jzg — a;ol « 1.37 fm. 

Our calculated form factors are shown in Fig. 3. 
is consistent with zero at = 0 and positive 
for all other values of . It rises with until it reaches 
a maximum value around 0.003, somewhere between 0.2 
and 1.0 GeV^, above which the data hint at a decrease. 
In the same range, is negative, with a decrease in 

magnitude with up to around 0.5 GeV^. Above that, 
the data are consistent with a constant. Note that data 
with the same spatial momentum transfer q = p — p^ 
tend to have strongly correlated errors, meaning that, 
e.g., the points below = 0-4 GeV^ form three clus¬ 
ters of correlated data that should not be interpreted as 
“bumpy” behavior in the form factor. The disconnected 
light-quark form factors have similar dependence on 
but their magnitude is two to three times that of the 
strange-quark form factors. 

We fit the Q^-dependence of each form factor using the 
z-expansion [27, 28], 


Ic 

'^max 

G(g2) = ^ akz\ 
k 


\/^cut V^^cut 

cut + + ^/tcut 


( 5 ) 


which conformally maps the complex domain of analytic- 
ity in to \z\ < 1. Although for physical quark masses 
the isoscalar threshold is fcut = at our pion mass 

we expect that the ui resonance is a stable particle be¬ 
low threshold; therefore, we use the isovector threshold 
Gut = (2TO.n.)^ in our fits. The intercept and slope of 
the form factor at = 0 can be obtained from the first 
two coefficients, ag (which we fix to zero for G^) and Oi. 
We impose Gaussian priors on the remaining coefficients, 
centered at zero with width equal to 5max{jaol, [oi]}. 
We truncate the series with k^ax = 5, but we have ver¬ 
ified that using fcmax = 10 produces identical fit results 
in our probed range of Q^. The resulting fit curves are 
shown in Fig. 3. 

The uncertainties are estimated as follows: 


1. Statistical uncertainties are computed using jack¬ 
knife resampling, using samples binned from four 
gauge configurations. These are shown as the inner 
error bands in Fig. 3. 

2. Fitting uncertainties are estimated by doubling and 
halving the widths of the fit priors, as well as per¬ 
forming the fits using different estimations of cor¬ 
relations in the data. 

3. For the light-quark disconnected form factors, we 
estimate excited-state errors by computing the 
form factors using a larger source-sink separation, 
1.37 fm. We assign the same relative error due to 
excited states to the strange-quark form factors. 

4. Uncertainties due to discretization effects are esti¬ 
mated by comparing against the form factors com¬ 
puted using the unimproved current. 

The outer error bands in Fig. 3 show the sum of 
these uncertainties in quadrature. We neglect finite- 
volume effects since they are highly suppressed by 
with = 5.9. Our resulting strange radii 

{r'^E mY — ~6dG|; j^^/d(3^lQ2^o and strange magnetic 
moment = G|,^(0) at pion mass 317 MeV are the 
following: 


(r|)" = -0.0054(9)(6)(11)(2) W, 

{rltY = -0.0147(61)(28)(34)(5) fm^, (6) 
= -0.0184(45)(12)(32)(1) 

where the four uncertainties are given in the same order 
as listed above, and is the nuclear magneton using 
the lattice nucleon mass, 1067(8) MeV. 

Although a controlled extrapolation to the physical 
point would require several lattice ensembles with vary¬ 
ing quark masses, an estimate can be made by com¬ 
bining strange-quark data with the equivalent obtained 
from the disconnected light-quark form factors. By it¬ 
self, the latter can be understood in the framework of 
partially quenched QGD [35] by introducing a third light 
quark flavor (which couples to the current in the quark- 
disconnected loop) and a bosonic ghost quark (which 
cancels all other loops). The dependence of the strange 
radii and magnetic moment on quark masses has been 
studied in SU(3) heavy-baryon chiral perturbation the¬ 
ory (GhPT) [36-38], and its partially-quenched general¬ 
ization [39-41]. At leading one-loop order, these observ¬ 
ables depend only on the mass mioop of a pseudoscalar 
meson composed of a nucleon valence quark and a quark 
from the vector current. For strange-quark and discon¬ 
nected light-quark observables, mioop is mK and m,r, re¬ 
spectively, and we can interpolate to the physical kaon 
mass. However, using typical values of the meson decay 
constant and meson-baryon couplings from phenomenol¬ 
ogy predicts a much stronger dependence on mioop than 
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FIG. 3. Strange-quark and disconnected light-quark electric and magnetic form factors, with statistical error bars. The curves 
result from the ^-expansion fits; the inner bands show the statistical uncertainty and the outer error bands show the combined 
statistical and systematic uncertainties (added in quadrature). The charge factors are not included. 
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FIG. 4. Linear combination of form factors, G% + r)G\f, 
probed by forward-angle parity-violating elastic ep scatter¬ 
ing experiments [5-7, 9-11, 13, 14]. The coefficient p depends 
on the scattering angle and ; for the lattice data we use the 
approximation p = AQ^, A = 0.94 GeV“^ [10]. In the low 
region we also show the linear dependence on resulting 
from the estimated charge radius and magnetic moment at 
the physical point. 

we observe, suggesting that the quark masses are too 
large for ChPT at this order. Therefore, we resort to a 
simple linear interpolation in We also adjust to 

the physical nuclear magneton, and obtain at the physi¬ 
cal point: 

(r%y = -0.0067(10)(17)(15) W, 

(^m)' = -0.018(6)(5)(5) fm^, (7) 

=-0.022(4)(4)(6) 


FIG. 5. Determinations of the strange magnetic moment: 
from direct lattice QCD calculations (this work and Ref. ]17]; 
red circles), models and phenomenology [16, 29-33] (green 
squares), and from a recent global analysis of parity-violating 
elastic scattering data [34] (blue diamond). 

where the first two uncertainties are statistical and sys¬ 
tematic (as estimated above). The third error is the dif¬ 
ference between the value at the physical point and on 
our lattice ensemble (using the physical nuclear magne¬ 
ton), and serves as an estimate of the uncertainty due to 
extrapolation to the physical point. 

The experiments run at forward scattering angles were 
sensitive to a particular linear combination of form fac¬ 
tors, G% + which we show in Fig. 4. Our results 

and the experimental data are both broadly consistent 
with zero, although the lattice data have much smaller 
uncertainties. This suggests that it will be quite chal- 
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lenging for future experiments to obtain a clear nonzero 
strange quark signal at forward angles. 

Figure 5 shows a comparison with some other deter¬ 
minations of the strange magnetic moment. The value 
from experiment has the largest uncertainty and is con¬ 
sistent with the other shown results. This work is in 
agreement with the other values within 2(T, except for 
two of the dispersion-theory scenarios [33], and has the 
smallest uncertainty. 

The techniques used in this work have proven effective 
in dealing with the longstanding problem of noise in dis¬ 
connected contributions to matrix elements. Although 
future calculations at near-physical quark masses will be 
needed to confirm our physical-point estimates, we have 
found very small contributions from strange quarks to 
proton electromagnetic observables: including the charge 
factor of —1/3 yields a roughly 0.3% effect in the proton 
fJ'r'i,, and 
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